Spin induced nonlinearities in the electron MHD regime 



Martin Stefan, Gert Brodin, and Mattias Marklund 

Department of Physics, Umea University, SE-901 87 Umed, Sweden 

We consider the influence of the electron spin on the nonlinear propagation of whistler waves. For this purpose 
a recently developed electron two-fluid model, where the spin up- and down populations are treated as different 
fluids, is adapted to the electron MHD regime. We then derive a nonlinear Schrodinger equation for whistler 
waves, and compare the coefficients of nonlinearity with and without spin effects. The relative importance of 
spin effects depend on the plasma density and temperature as well as the external magnetic field strength and 
the wave frequency. The significance of our results to various plasmas are discussed. 
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I. INTRODUCTION 

The dynamics of dense ionized matter, cold matter in strong 
magnetic fields, and nonlinear degenerate plasmas has appli- 
cations to both laboratory and naturally occurring systems. 
Many such matter states goes under the collective notation of 
quantum plasmas. A multitude of studies devoted to quantum 
plasma effects can be found in the literature, much of it in- 
spired by works like e.g. Refs. UI2]. In the last decade there 
has been somewhat of a sur ge in the interest of quantum plas- 
mas d, 0, H, H, 0, H H, [roTfTll [l2ll . Interesting applications 
of this field can be found in for example plasmonics fl3l[l4ll . 
quantum wells lfl5tl and ultracold plasmas iflrill . Common to 
such applications are rather "extreme" parameters, compared 
to most laboratory and space plasmas. More specifically, the 
plasma densities needs to be very high and/or the temperatures 
correspondingly low. For astrophysical plasmas, the situation 
is somewhat different since the strong magnetic fields lfl7i[l8ll 
may induce various types of quantum effects. In most of the 
above studies, the spin effects plays little or n o dy namic role. 
The inclusion of collective spin dynamics 20]gives rise to 
new modes in plasma, both at fluid ||9l |2l|] and kinetic scale 
fl2tl . Indeed, even dusty plasmas can show interesting mag- 
netization effects JH, In Ref. |Q~D] the picture outlined 
above, concerning the necessary parameter space for quantum 
effects to be important, was to some extent modified, as it 
was shown that the spin properties of electrons can be impor- 
tant in plasmas even outside the high density/low temperature 
regime, also for moderate magnetic field strengths. Moreover, 
a recent focus on the nonlinear regime in quantum plasmas 
ll24l I25I l26ll makes the question of magnetization nonlineari- 
ties interesting. 

Motivated by the above, we will in the present work further 
extend the analysis put forward in Ref. Ill ill . In that work, the 
electrons were described using a two-fluid model, where the 
spin-up and spin-down populations relative to the magnetic 
field were treated as different fluids in the standard magneto- 
hydrodynamic (MHD) regime. Here we will extend that treat- 
ment to cover the electron-MHD (EMHD) regime. In particu- 
lar we will study weakly nonlinear whistler waves, and derive 
a nonlinear Schrodinger (NLS) equation for the slowly vary- 
ing amplitude, both using classical weakly relativistic theory 
and the recent two-fluid spin model, adapted for the EMHD 
regime. By comparing the nonlinear coefficients in the dif- 



ferent models and their dependence on the plasma parame- 
ters (temperature, density external magnetic field strength), 
the relative importance of the electron spin effects in vari- 
ous regimes can be deduced. The result that electron spin ef- 
fects can be important in other regimes, as compared to certain 
much studied quantum effects [such as the Bohm-de Broglie 
potential (see also Ref. l27ll for a discussion) and the Fermi 
pressure], is confirmed. Finally we compare the relative im- 
portance of electron spin effects in the EMHD regime with 
that in the standard MHD regime. 



II. TWO-FLUID MODEL WITH SPIN 

The purpose of our work is to compare the nonlineari- 
ties from classical and quantum effects, respectivel y, in the 
EMHD regime. As a starting point, we follow Ref. 12811 and 
derive the EMHD model. We assume that the time-scale of 
interest is small enough that the ion motion, due to their large 
mass compared to the electrons, can be neglected. We also ne- 
glect the displacement current in Amperes law. Furthermore 
by assuming an isothermal pressure model and no dissipation, 
the governing equation can be written as 

(B — dgAB) = -aVx [(V x B) x (B-^AB)] , (1) 

where B is the magnetic field, d 2 — c 2 / (0 2 e , (0 2 e — q 2 n/m e £o 
being the electron plasma frequency, and a = l/nq e jj,Q. 
Studying the linear modes of this equation propagating paral- 
lel to the magnetic field one finds whistler waves with disper- 
sion relation co(k) = (O ce c 2 k 2 / 0) 2 e . However, since this model 
does not allow any density fluctuations, an attempt to derive an 
NLS-equation shows that the model do not give rise to any cu- 
bic nonlinearities. To still be able to do the intended compari- 
son, we still use the assumptions corresponding to the EMHD 
regime but perform a more general treatment, allowing for rel- 
ativistic particle velocities and using a multifluid model that 
permits density perturbations. 

The quantum model used in the comparison is obtained for- 
mally by starting from the Pauli Hamiltonian as is done in Ref. 
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lID], using ensemble averaging to obtain the fluid equations 
d 



dt 
d 



n + V-(nv) = 



(2) 



m„n\ — +v-V)v = ^n(E + vxB)-Vp + F spin (3) 



where 



F spin = ±AiB«VB, 



(4) 



the sign depending on the spin orientation S± = ±(n/z/2)B 
relative to the magnetic field. Thus, assuming a two-electron 
fluid model, for which the distinction between the fluids is 
through their relative spin orientation, we have 

d\± 

m e n± = n±q e (E + v± X B) ± jl B n±'VB — k B TVn±(5) 

dt 



dt 



-V-(n±v±), 



(6) 



where the subscript ± denotes spin orientation parallel or 
anti-parallel to the external magnetic field respectively, and 
B = |B|. In case the up- and down spin populations are not 
equal, there will be a net magnetization and a corresponding 
magnetization current. Thus within this model the total cur- 
rent density to be used in Ampere's law is written 



j = q e (n + y+ +n_v_) + /x^V x (n+B — n_B) 



(7) 



where the last term is the magnetization current due to the 
spin, and B = B/B is a unit vector in the direction of the mag- 
netic field. 



HI. LINEAR THEORY 

Linearizing the momentum and fluid equations around a 
constant magnetic field Bo = Bqz, and assuming transversal 
waves propagating parallel to this external magnetic field we 
can deduce 



from Eq. (0, where 



<7.v 



AT 



v = ctE, 



m s {(0 2 s - (o 2 ) 



(8) 



(Or, 



a> C s 

-i(0 



is the conductivity tensor, and <o cs = q s Bo /m s is the cyclotron 
frequency for particle species s, and we let the vectors here 
just contain the parts perpendicular to z. Note that since 
the variations of B is nonlinear in the amplitude for parallel 
propagation, the spin effects do not enter here. Furthermore, 
when the thermal energy k B T is much larger than the energy 
difference between the spin states, /i b Bq, the difference be- 
tween the number density of the spin-up and down popula- 
tions, «q + — «o_, in the thermodynamic ground state is expo- 
nentially small (proportional to exp(— fl B Bo /k B T)), and hence 



we can omit the linearized part of the magnetization current. 
Thus in this approximation no quantum effects remains in the 
linearized theory. From Maxwell equations, we then obtain 



DE = 



(9) 



were 



D = 



k 2 c 2 , . v , 

(0 2 ^ 



03 



ps 



-CO —l(O cs 
'f (0 ((0 2 - CO 2 ) V i( °cs -03 



Thus, we obtain the general dispersion relation 



oo 2 = k 2 c 2 



to; 



1- 



(Ocs 
(O 



-1 



(10) 



(11) 



for this geometry 

To obtain the EMHD limit, we disregard the displacement 
current in Amperes law and regard the ions as fixed. This 
corresponds to neglecting the unit matrix in Eq. (O and letting 
(O ce *C (0. By this procedure we get the dispersion relation 



(0(k) 



^c 2 k 2 



Bpk 2 
q e n e lk> 



= aB Q k z 



(12) 



for small amplitude whistler waves propagating parallel to the 
external magnetic field. 

The dispersion relation <fT~2b was obtained by starting from 
the two-fluid model and then taking the EMHD limit. If one 
would start directly from the EMHD plasma equation, the cor- 
responding result would be 



(o(k) 



aB Q k 2 



1 



: 2 k 2 /(0 2 ' 



(13) 



This difference is due to the fact that for the EMHD assump- 
tions to apply, the parallel (to the external magnetic field) 
wavenumber k\\ must obey k\\ -C (O p /c. Thus in our case with 
parallel propagation (k =kn) we must approximate the de- 
nominator with unity, in which case the different expressions 
agree. As a side note, for general directions of propagation 

the factor ( 1 + c 2 k 2 / (O 2 ) appears correctly from EMHD the- 
ory, but since only the perpendicular part of the wavenumber 
k± is allowed to be comparable to the inverse skin depth (i.e. 

k ± ~ <o p /c), we have (1 +c 2 k 2 /<o 2 y = (1 +c 2 k\/(0 2 y in 
that case. 

As usual, by using an ansatz of a weakly modulated ampli- 
tude in the EMHD model Eq. (Q~|), neglecting nonlinear terms 
and higher order dispersion we obtain the linear part 







v' g d 2 - 






+ 2 dz 2 _ 



B x = 0, 



(14) 



of the one-dimensional NLS equation. Here v g is the group ve- 
locity and v' = dv g /dk is the group velocity dispersion. Since 
the spin effects do not enter linear theory, these coefficients 
are the same in our classical and quantum mechanical mod- 
els. 
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IV. CLASSICAL NONLINEAR THEORY 

To explore nonlinearities due to relativistic effects, the mo- 
mentum equation is modified by letting \ s — > yv v in the left 
hand side, and it thus reads 

n * ( T + v * ' V ) rr* = — n * ( E + v * x B ) - v ?t v "^ ( 15 > 
V ot J m s 

where y = 1/(1 — v 2 /c 2 ) and we have introduced the thermal 
velocity v st = (kgT /m s ) 1 ^ 2 for species s. The y-factor can be 
Taylor expanded to first order and will thus result in purely 
cubic nonlinearity in the velocity. Including this nonlinearity 
only, it is straightforward to deduce the NLS equation 



d_ 
Tt 



2a 3 B Q k 6 



s dz 2 dz 2 



co 2 I 



-c 2 k 2 /co 2 ) 



\B X \ Z B x = 0. 
(16) 

The nonlinear term above will be complemented by nonlin- 
ear density modifications induced by the ponderomotive force. 
Within a model that only includes the electron dynamics, the 
density modifications will be limited due to the general ten- 
dency of charge neutrality. However, for sufficiently long 
pulses the low frequency ion dynamics will start to contribute 
to the nonlinear behavior of the electrons, and it turns out that 
a fair comparison between quantum and classical nonlinear- 
ities must include this effect. To capture the ponderomotive 
nonlinearities, we start with Eq. (fl3T > (for simplicity omitting 
the relativistic contribution, that we already know), for a two- 
fluid ion-electron model . Again neglecting the displacement 
current in Amperes law, linearizing as previously and using 
the Maxwell equations we obtain the system 



det [£>, 



op J 



/ VeK \ 
E\f 
"elf 
Vflf 

V «»if / 



I name dz 1 



Adj [A 



op J 



\B X \ 2 \ 






-2 d 



\ 







\B, 



(17) 



for the low-frequency variables, where 



£> op = 










rile 
d 
dz 





m e ~d~z 

i 










"4 



\ 

£o 


kJi d 
mi ~d~z 
d 

Tt 



(18) 



and we can read off that 

Vflf 



= kJB 



= Ki\B x 



2 



(19) 
(20) 



The coefficients K e and jc, are determined by solving the corre- 
sponding differential equation, obtained from Eq. (fTTJi, using 
Greens function techniques, and the result is: 

(0 2 p Jm e + (0 2 pe / mi 



-2v 



noHo 



aUvl-vfJ + aUvl-v 2 ) 



(21) 



Now that the low frequency perturbations have been deter- 
mined, the back reaction on the original time scale can be cal- 
culated. We note that on this fast time scale we return to ne- 
glecting ion motion. Then, we obtain 



V x 



noveif 



d ( 1 



dz \Hoq e n d 
k 3 k 2 B, 



zxB 



n e \ f \ x B 



JO 



(22) 



kJBJ-B^p 



where the subscript x indicates the i-component of the vector. 
Combining this result with the linear theory we obtain an NLS 
equation that reads 



d d . 

Tt +v 'Tz ]£ - 



v' d 2 B x 



2 dz 2 
q e B()k 2 c 2 /m e v g — k^c 2 



<(1 



-c 2 k 2 /co 2 



pel 



(23) 



\B X \% 



V. FULLY NONLINEAR THEORY 

As opposed to the case of a normal EMHD plasma, in the 
quantum case we have one equation for each electron spin 
direction, and the extra term ±/Xb«±VB due to the spins in- 
fluence on the magnetization. One can note that if the spin 
populations are exactly equal in density, when adding the two 
force equations this term will vanish, and this corresponds to 
the classical case. However, if there is a slight difference in 
number density, nonlinear fluctuations in the magnetic field 
will be induced. To explore this effect, we try to derive an 
EMHD model , but now using the two-fluid spin model. Eq. 
(Q]i is then replaced by 



2 Nq e ^ = ^ V x [(V x B) x B] - 2/i B V x [(V x (nB)) x B] 



-2pi B V x (nVB) , 



(24) 



where the average number of electrons = («+ + n_)/2 tend 
to be deviate little from the unperturbed density (due to charge 
neutrality), but the difference between the electron species n = 
(n + — n_)/2 can vary more. Furthermore we introduce the 
notation V = v+ + v_ and v = v + — v_. We here point out 
that in the approximation considered, where the unperturbed 
density difference difference is neglected, similar as before 
the linear treatment give agreement with the classical case. 

Next, to calculate the low frequency perturbations of n, we 
need the difference between Eqs. © for the two species, 
which is written 



0= -[(Nv + nV) xB] 



-2N 



B 



B 



2kgT dn 



(25) 



q e B dz x yJ q e dz ' 
Filtering out the low frequency time scale we obtain: 



at, 



Mf = - 



1 



fiomv 2 . dz 2 



Br 



IB* 



(26) 
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and 



n xf = 2N 



k B TB 



\B x \ Z + \By\ 



(27) 



where we have introduced the electron thermal velocity, v te = 
(ksT /m e ) l l 2 . Due to the reduced geometry of the problem, 
with parallel propagating circularly polarized modes, all sec- 
ond harmonic density perturbations can be neglected, and all 
second order nonlinearities in the magnetic field also vanish. 
Furthermore, N\f <C «if as a consequence of the system tend- 
ing towards charge neutrality, and thus only «if needs to be 
considered here. 

Inserting this now in Eq. (f24b and considering the original 
time scale the only nonvanishing contribution to the nonlinear 
constant is 

- MbV x [(V x (nft)) x B] = !6iN-^- \B X \ 2 B X . (28) 

Thus, the full NLS equation including all the effects discussed 
above (relativistic nonlinearity, classical density perturbations 
induced by the ponderomotive force, and a spin dependent 
density modification, driven by the nonlinear magnetic dipole 
force) will be (cf. Eq. £23)) 



d d , 

2a^B k 6 



v' g d 2 B x | 4^i 2 k 2 
2 dz 2 kgTBoqt 



\B X \ Z B X 



-^B k 2 c 2 - k 3 c 2 



\B x \ 2 B x = 0. 
(29) 



This equation is the main result of this paper. From the mag- 
nitude of the nonlinear coefficient, one can determine the 
regimes in which the spin terms can dominate and be respon- 
sible for e.g. soliton formation. 



VI. DISCUSSION 

In the present paper we have studied weakly nonlinear 
whistler waves propagating along the magnetic field. A non- 
linear Schrodinger equation has been derived for the case of 
classical nonlinerities (see Eq. (|29l). The nonlinear coef- 
ficient than gets two contributions; from relativistic effects 
and from low-frequency density modifications induced by the 
ponderomotive force. Taking spin effects into account, within 
in a electron two-fluid spin model, it is found that the low- 
frequency part of the magnetic dipole force separates the spin 
up and spin down populations. Due to the different magneti- 
zation currents from the two populations, a spin contribution 



to the nonlinear coefficient then arises. Firstly, comparing the 
spin contribution to the nonlinear coefficient with the relativis- 
tic contribution, we see that the former is larger provided that 



miti 2 (0j e 



k B Tm 2 C\ 



(30) 



Here we have used the lowest frequency allowed by the model 
CO ~ CO C i to get a condition that is relatively easy to fulfill. 
However, we must also compare the spin induced nonlinear- 
ity against the contribution from the nonlinear density oscilla- 
tions induced by the ponderomotive force. It is found that the 
former dominates when 



1 < 



TlCOce TlCOc, 



•iC 2 



(31) 



where we have used that the maximum value of kc is roughly 
(Ope, due to the limitations imposed by the geometry in 
combination with the EMHD approximation. The factor 
h(O ce /miC^is the condition for nonlinear spin effects to domi- 
nate, when a similar comparison is made in the standard MHD 
regime, according to Ref. [11]]. To get a more favorable com- 
parison (i.e. a condition that is easier to reach under laboratory 
conditions) than in these previous works, the second factor, 
hcOce/mjvjj, must be larger than unity. This is unfortunately 
not the case for the parameters usually found in laboratory 
conditions. However, astrophysical plasmas with parameters 
fulfilling both conditions (130b and OTb can be found, e.g. in 
the vicinity of pulsars or magnetars [17], and thus we note that 
effects associated with the electron spin can be more impor- 
tant than the classical relativistic and ponderomotive nonlin- 
earities in such environments. 

The present study has focused on the EMHD regime. While 
it is shown that spin effect certainly can be important during 
e.g. astrophysical plasma conditions, our study suggest the 
standard MHD regime 111 111 can be more affected by the elec- 
tron spin properties during laboratory conditions . However, 
much more work remains to be done in order for this con- 
clusion to be settled, as the picture may change when a more 
general geometry is considered, or when kinetic effects lfl2ll 
are taken into account. 
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